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Abstract 

Charge separation is an important consequence of the Chiral Magnetic 
Effect. Within the framework of a multi-phase transport model, the ef- 
fects of final state interactions on initial charge separation are studied. We 
demonstrate that charge separation can be significantly reduced by the evolu- 
tion of the Quark-Gluon Plasma produced in relativistic heavy ion collisions. 
Hadronization and resonance decay can also affect charge separation. More- 
over, our results show that the Chiral Magnetic Effect leads to the modifica- 
tion of the relation between the charge azimuthal correlation and the elliptic 
flow that is expected from transverse momentum conservation only. The 
transverse momentum and pseudorapidity dependences of, and the effects of 
background on the charge azimuthal correlation are also discussed. 
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1. Introduction 

Charge separation along the angular momentum direction has been inves- 
tigated in relativistic heavy ion collisions [H, lij- The experimental study was 
motivated by the theoretical investigation of the Chiral Magnetic Effect 0- 



lOj . The Chiral Magnetic Effect is related to the fact that the hot and dense 
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matter created in heavy ion collisions can form V and CV odd metastable 
domains where the parity and time-reversal symmetries are locally violated. 
In the early stage of a non-central relativistic heavy ion collision, the mag- 
netic field can reach a magnitude on the order of 10^^ T. In the presence of 
a strong magnetic field, these topologically non-trivial domains impose con- 
straints on quark chiralities and induce a separation of negative and positive 
particles in the direction of magnetic field (i.e. system angular momentum). 
In spite of large theory uncertainties, the experimental results are consistent 
with Chiral Magnetic Effect expectations. 

In addition to the Chiral Magnetic Effect, other effects can also con- 
tribute to charge separation and/or charge correlation. Bzdak et al. found 
that the contribution due to transverse momentum conservation is compara- 
ble in magnitude to the prediction of the Chiral Magnetic Effect as well as 
the data Wang also argued that the measured data can be accounted 
for by cluster particle correlations and new physics may not be required to 



explain the data [12[. Schlichting and Pratt argued that local charge con- 



servation, when combined with elliptic flow, explains much of experimental 



measurements [13[. To our knowledge, no previous studies have included 
the dynamical effects of final state interactions, such as parton cascade and 
resonance decay, on the experimental charge separation observable. On the 
other hand, these final state interaction effects have been found important 
for many experimental observables, such as elliptic flow and particle yields. 
In the following, we will address the problem of whether an initial charge 
separation will be able to survive the final state interactions. 

The paper is organized as follows. In Section [21 we give a brief descrip- 
tion of A Multi-Phase Transport (AMPT) model and its adaptation for the 
study of charge separation. Results on the charge separation observable are 
presented in Section E] followed by a summary in Section HJ 

2. A Multi-Phase Transport model 



The AMPT model |14|-|16[| is a dynamical transport model that includes 
four different stages in relativistic heavy ion collisions: the initial condition, 
partonic interactions, the conversion from partonic matter into hadronic mat- 
ter, and hadronic rescatterings. The initial condition, which includes the 
spatial and momentum distributions of minijet partons and soft string ex- 
citations, is obtained from the Heavy Ion Jet INteraction Generator (HI- 



JING) model 17|, [iSj. There are two options for doing the parton evolution 
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and hadronization in the AMPT model. One option is the default model 
which includes only interactions of minijet gluons via Zhang's Parton Cas- 
cade (ZPC) jioj and uses the Lund string fragmentation model 2^ to turn 
partons into hadrons. The other option is the string melting model. It starts 
the parton evolution with a quark-anti-quark plasma from the dissociation 
of strings. It recombines partons via a simple coalescence model to produce 
hadrons 22|. Dynamics of the subsequent hadronic matter is then described 
by A Relativistic Transport (ART) model 2l|. The default model has good 



agreement of particle spectra with experimental data but it significantly un- 
derestimates the elliptic flow. In contrast, the string melting model can only 
describe low transverse momentum spectra, but the agreement with the el- 
liptic flow data is much better. In addition, the AMPT model has been used 
to study other observables, such as strangeness |23l. |24|] . charm 25|, J/^|_ pro- 



duction |26l-l28j . two-pion correlation function [29|, dijet correlations [30 
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triangular and higher order flows [32|, |33 



The AMPT model with string melting starts with a quark-anti-quark 
plasma and it will be used for the study of the effects of flnal state inter- 
actions on charge separation. The Chiral Magnetic Effect is not built into 
the AMPT model. In order to separate a fraction of the charges initially, 
we switch the Py values of a fraction of the downward moving u quarks with 
those of the upward moving u quarks, and likewise for d and d quarks. The 
coordinate system is set up so that the x-axis is in the reaction plane and 
the y-axis is perpendicular to the reaction plane with the 2;-axis being the 
incoming direction of one nucleus. The above procedure ensures total mo- 
mentum conservation while giving momentum kicks to produce an upward 
initial current. The current implementation of the ART model does not con- 
serve the electric charge. In the following, we turn off the hadron evolution. 
The strong parton cascade provides the most contribution to the evolution 
and the exclusion of hadron evolution in the string melting model is not ex- 
pected to make signiflcant changes to the flnal results. Resonance decays are 
implemented to ensure charge conservation and are included for the study of 
charge correlations. 



3. Charge separation in heavy ion collisions 

To measure charge separation possibly coming from local strong parity 
violation in relativistic heavy-ion collisions, the STAR experiment studied a 
charge azimuthal correlation observable (cos(0q + 4>/3 — as proposed 
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Figure 1: Centrality dependence of {cos{<j)a + </>^)) in Au+Au collisions at ^5^=200 GeV 
(with a 10 mb parton cross section). The different symbols represent different percentages 
of initial charge separation in AMPT calculations. The stars represent experimental data, 
where the two surrounding curves give the systematic uncertainty for data. Some points 
are slightly shifted for clarity. 

by Voloshin [3^]. Here, a and /3 represent the signs of electric charges and 
can be positive or negative, while "^rp is the azimuthal angle of the reaction 
plane. By measuring the correlations for same-charge and opposite-charge 
pairs, the data show some hints of charge separation, which is consistent with 
the expectation of the Chiral Magnetic Effect jl[ lij. In the following, we will 
look at how the initial charge separation contributes to the charge azimuthal 
correlation. The reaction plane azimuthal angle will be set to zero degrees. 
Au+Au collisions at ^/s^ = 200 GeV with a 10 mb parton cross will be 
studied and light quark and charged pion correlations will be analyzed. 

Fig. [1] presents the charge azimuthal correlation as a function of central- 
ity from the AMPT simulations. Since the initial charge separation could 
depend on centrality, different percentages of initial charge separation are 
used for each centrality bin to look for possible centrality dependence. The 
STAR correlation data are also shown for comparison. For the same-charge 
correlation, results from the AMPT model without initial charge separation 
have smaller magnitudes than data. As the percentage of initial charge sep- 
aration increases, the magnitude of the correlation increases. The increase is 
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Figure 2: Centrality dependence of (cos((/)q, + (/i/g)) for different stages in AMPT cal- 
culations with an initial charge separation percentage of 10% for Au+Au collisions at 
V'^NN=200 GeV (with a 10 mb parton cross section). Some points are slightly shifted for 
clarity. 

not linear in the initial charge separation percentage. A percentage of 10% 
for initial charge separation can describe well the STAR measurements. For 
the opposite-charge correlation, it seems that initial charge separation is not 
necessary for all centralities except the most peripheral bin of 60-70%. For 
the centrality bin of 60-70%, a percentage of 10% is indeed needed to match 
the experimental observation. In other words, the observed opposite-charge 
correlation changes much faster than the AMPT results. Even though 10% 
initial charge separation can describe both the same-charge and opposite- 
charge correlations for the 60%- 70% centrality bin, it is difficult to describe 
the centrality dependence of both the same-sign and opposite-sign correla- 
tions with initial charge separation alone. We also notice that the results 
with 5% initial charge separation are almost identical to those without ini- 
tial charge separation. This indicates that it might be challenging to observe 
an initial charge separation of 5% or less in the presense of strong final state 
interactions. 

To understand how the charge correlation observable evolves in heavy- 
ion collisions. Fig. [2] shows the centrality dependence of charge correlations 
for different stages in the AMPT model. With an initial percentage of 10% 
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charge separation, the initial charge correlations are quite large (solid and 
dash lines), with a magnitude of about 5.5 x 10~^. After strong parton 
cascade, charge correlations are significantly reduced especially for central 
collisions because of frequent parton interactions under high parton density. 
The charge correlations are recovered partly from hadronization as coales- 
cence reduces the number of particles while combining quarks into hadrons. 
Resonance decays act opposite to coalescence and reduce charge correlations 
in the hadronic phase. The final charged pion correlations have magnitudes 
comparable with those of final partons. Related to the charge correlation is 
the percentage of charge separation. Its centrality dependence has the same 
qualitative evolution where parton cascade and resonance decay decrease 
while coalescence increases the percentage. From a percentage of charge sep- 
aration of 10% in the beginning, only 1-2% percentage remains at the end 
with more peripheral collisions having larger percentages. 

As a comparison, the charge correlations at different stages with no initial 
charge separation are shown in Fig. |3]for different centrality bins. Before the 
parton cascade, both the same-charge and the opposite-charge correlations 
are consistent with zero. After the parton stage, both correlations become 
negative with the same-charge correlation having the larger magnitude. Neg- 
ative correlations indicate that the correlated charges move together and they 
are not separated. Coalescence increases the magnitude for the same-charge 
correlation and resonance decay decreases it as in the case with a non-zero 
initial charge separation. However, for opposite charges, coalescence reduces 
the correlation. If the opposite-charge correlation is calculated including 
charged rho mesons in addition to charged pions, it has a magnitude that 
is larger than that of quarks after parton cascade. This shows that when 
there is no initial charge separation, the opposite-charge correlation from 
coalescence is not equally distributed among different species combinations. 

Recently, Bzdak et al. found that transverse momentum conservation can 
contribute to the charge correlations with magnitudes comparable to exper- 



imentally observed correlations [ll[. The charge correlations can be calcu- 
lated from transverse momentum conservation alone. Under the assumption 
that all particles have the same average transverse momentum, there is a 
simple relation between the charge correlations and the elliptic flow. Both 
the same-charge and opposite-charge correlations are equal to -V2/N for suf- 
ficiently large N. Here V2 is the elliptic flow coefficient and N is the total 
number of produced particles (similar results were also obtained in jssi 36|). 
The opposite-charge correlation can be affected by factors other than trans- 
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Figure 3: Centrality dependence of {cos{(pa + 4>p)) for different stages from AMPT calcu- 
lations without initial charge separation for Au+Au collisions at ^^5^=200 GeV (with a 
10 mb parton cross section). Some points are slightly shifted for clarity. 

verse momentum conservation. In the following, we will look at how this 
relation holds for the same-charge correlation. Fig. H] presents the same- 
charge correlation as a function of —V2/N for different stages from AMPT 
calculations without (0%, open symbols) and with (10%, solid symbols) ini- 
tial charge separation. Here f 2 and are for particles with pseudorapidity 
\ri\ <1. The AMPT results without initial charge separation are consistent 
with the expectation that (cos(0q, + ^/j)) = —V2/N, which is shown in the 
figure by a dashed line. It indicates that the same-charge correlation is driven 
by transverse momentum conservation in the AMPT model without initial 
charge separation. On the other hand, the AMPT results with initial charge 
separation are much lower than the expected relation. It is interesting to see 
that the linear relation between the same-charge correlation and —V2/N is ap- 
proximately preserved with a coefficient much large than 1. Since the AMPT 
results with 10% initial charge separation can describe the same-charge cor- 
relation data well as shown in Fig. [H transverse momentum conservation can 
only partly account for the measured charge correlation data. 

In more detail. Fig. |5] shows the dependences of charge correlations on 
the average of the transverse momentum (p+ = {pt,a + Pt,p)/'^) of two final 
charged pions for the 30-50% centrality bin. For the same-charge correla- 
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Figure 4: (cos(0q + (f)^)) as a function of —V2/N for different stages in AMPT calculations 
without (0%, open symbols) and with (10%, solid symbols) initial charge separation for 
Au-|-Au collisions at ^s„„=200 GeV (with a 10 mb parton cross section). The dashed 
line represents the relation of (cos((/)q + ^/s)) = ~V2/N . 



tion, the magnitudes of results from the AMPT model without initial charge 
separation are smaller than those of data while 10% initial charge separation 
can increase the magnitudes to reproduce data. For opposite-charge pairs, 
the correlation with no initial charge separation is consistent with data while 
the correlation with 10% initial charge separation increases weakly with p+ 
and is a little higher than data. All of these results are consistent with the 
integrated correlations which are presented in Fig. [1] From transverse mo- 
mentum conservation, (cos((^o + 0/?)) is proportional to with n = 2 to 



3 ll| . The curves in Fig. [5] show the power- law fits to the same-charge cor- 
relations from the AMPT model. The power n is 2.24 ± 0.27 when there 
is no initial charge separation, consistent with the expectation from trans- 
verse momentum conservation. However, the power n deceases to 1.54 ± 
0.18 when 10% initial charge separation is included. In addition, we found 
that the charge correlations depend very weakly on p_ = \pt,a — Pt,/?!- In 
particular, the opposite-charge correlation increases gradually to a level of 
about 0.1 X 10~^ while the same-charge correlation stays at a constant level 
of about 0.25 x 10~^. Even though the magnitude of the same-charge correla- 
tion is smaller than the experimental data, the integrated value is consistent 
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FigureS: (cos(0q +0^)) as a function of p+ = (pt,a+pf^^)/2 in AMPT calculations without 
(0%, diamonds) and with (10%, squares) initial charge separation for the 30-50% centrality 
bin in Au+Au collisions at ^s„„=200 GeV (with a 10 mb parton cross section). The curves 
are power-law fits for the same-charge correlations from the AMPT calculations without 
(dash) and with (solid) initial charge separation, and the stars represent experimental 
data. Some points are slightly shifted for clarity. 

with experimental data because the lowest p_ bin carries the highest weight. 

Fig. [6] presents charge correlations as functions of the pseudorapidity dif- 
ference {At] = \ria — rii3\) of two final charged pions for the 30-50% centrality 
bin. Again we see that results from the AMPT model without initial charge 
separation can describe the opposite-charge correlation, while initial charge 
separation is needed to reproduce the same-charge correlation data. It is 
worth mentioning that the strong dependence on the pseudorapidity differ- 
ence cannot be obtained in the present calculations from transverse momen- 
tum conservation [ll|. More realistic longitudinal dynamics in the AMPT 
model contributes to the better description of the dependence on the pseu- 
dorapidity difference. 

In addition to the angular correlation (cos(0a + (p/s)), charge separation 
also shows up in the angular correlation (cos(0a — 0^))- The former is free 
of reaction plane independent backgrounds while the latter is also sensitive 
to reaction plane independent backgrounds. Charge separation increases the 
opposite-charge correlation and decreases the same-charge correlation for the 
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Figure 6: (cos((/)q + 0^)) as a function of Ary = \ria — Vpl in AMPT calculations without 
(0%, diamonds) and with (10%, squares) initial charge separation for the 30-50% centrality 
bin in Au+Au collisions at ^s„„=200 GeV (with a 10 mb parton cross section). The stars 
represent experimental data. Some points are slightly shifted for clarity. 

former, while it decreases the opposite-charge correlation and increases the 
same-charge correlation for the latter. We will look at the centrality de- 
pendence of the charge correlation (cos(0q, — ^/j)) in Fig. [7l When there is 
no initial charge separation, the AMPT results have the same trends as the 
experimental data for both the same-charge and the opposite-charge corre- 
lations. However, the correlations are much lower than those observed ex- 
perimentally. Charge separation brings the same-charge correlation closer to 
data and the opposite-charge correlation farther away from data by amounts 
comparable to those for (cos(0q, + ^/j)). However, the changes are not enough 
to make up for the large difference between the AMPT results and the ex- 
perimental data. Additional backgrounds that can significantly increase the 
correlations are needed in order to describe the data. 

4. Conclusions 

In summary, final state interactions play an important role on charge 
separation in relativistic heavy-ion collisions. Parton cascade and resonance 
decay significantly reduce the charge separation from 10% in the initial state 
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Figure 7: Centrality dependence of (cos((^a — 0^)) from AMPT calculations without 
(0%, diamonds) and with (10%, squares) initial charge separation in Au+Au collisions 
at y^s^=200 GeV (with a 10 mb parton cross section). The stars represent experimental 
data. Some points are slightly shifted for clarity. 



to 1-2% in the final state. Therefore, it is essential to take these final state 
effects into account for studies related to charge separation. Our results also 
suggest that mechanisms beyond transverse momentum conservation will be 
needed even for the description of the same-charge correlation. 

Our approach includes the effects of local charge conservation and trans- 
verse momentum conservation automatically. However, detailed magnetic 
field evolution [37|, or fluctuating domain sizes, or different topological charges 
are not included. These effects can lead to different charge separation per- 
centages for different centralities. But they are not likely to help improve 
the simultaneous description of both the same-charge and opposite-charge 
correlations, and both (cos(0q, + 0/3)) and (cos(0a — 0^)). Schlichting and 
Pratt recently demonstrated that charge balancing can affect the difference 



between the opposite-charge and same-charge correlations [13|]. This and 
other possible mechanisms certainly deserve further study for a satisfactory 
understanding of experimental data. 
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